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Abstract 



We extend a modal theory of diffraction by a set of parallel fibers to deal 
with the case of a hard boundary: that is a structure made for instance of 
air-holes inside a dielectric matrix. Numerical examples are given concerning 
some resonant phenomena. 

I. INTRODUCTION 

Photonic crystals (PCs) are periodically modulated structures that present the properties 
of having photonic band gaps In the case of 2D photonic crystal where it is possible 

to study separately s- and p-polarized fields, it is now known that large gaps are easier to 
obtain with inverted-contrast crystals (i.e. air holes in a dielectric matrix) with p-polarized 
fields than with s-polarized fields In this work we extend a multi-scattering theory by 
cylinders in order to deal with inverted-contrast crystals with a hard boundary, that 

is a structure made of dielectric or metallic inclusions in a dielectric matrix embedded in 
vacuum (fig. 1). This theory allows to quantify the importance of the hard boundary when 
the device is considered as a 2D PC and also to study the propagation phenomena when the 
device modelizes a photonic crystal fiber [|T0| -[T3||. 



II. THEORY OF DIFFRACTION 

We start by constructing the generalized scattering matrix of a set of N dielectric or 
metallic rods embedded in a dielectric cylinder of circular cross section VL (fig.l) and ra- 
dius R. We use a Cartesian coordinate system {0;x,y,z) with origin at the center of fl. 
The electromagnetic fields considered here are harmonic fields with a time dependence of 
exp{—iut). 

The rods are denoted by {-D„}, they are filled with a dielectric (relative permittivity En) 
or metallic material. The cylinder in which the fibers are contained has relative permittivity 
Er, finally this cylinder is embedded in vacuum. Due to the invariance of the medium along 
the 2;-direction we look for solutions of Maxwell equations with a z dependence of exp {i'yz) . 
Under this assumption, it is easily shown that all the components of both magnetic and 
electric fields are known once and Hz are known. Denoting 
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the following propagation equation is satisfied 

A^F+x'F =0 

we denote Aj_ = ^ + = {^) ~ 7^ and k is the wavenumber. We denote 

Xl = k'^Sr - 7^ 
Xl = k'^Sn - 7^ 
xl = k^eo - 7^ 

The total exterior field F is expanded in Fourier-Bessel expansion outside Q 
F (r, ^) = J2 K^^- (Xor) + Ft^Hi^^ (xor)] e^'"^ |r| > R 

m 

Remark: J2m '^m {xo^) e*™''^ represents the incident field. 

The total field F inside fl writes in the vicinity of the boundary dQ of Q 



(1) 



F (r, ^) = J2 iXrr) + F^fHi^^ ^^^r)] e 



irrup 
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We denote 

so that the transmission conditions on dVt write 

_|_ ^0!.+ = - _|_ 



_|_ pti,- 



+ F° 
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where are boundary impedance operators easily deduced from Maxwell equations. Con- 
ditions flSD lead to 




S2 
Si S2 
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where are linear operators deduced from (|^) linking the incoming and outgoing parts of 
the fields on the boundary of Vt. 

Remark: In case where the boundary is not circular the theory can still be applied 
provided the expansions are restricted respectively to the smallest circle of center O 
containing Vt and to the greatest circle contained in Vt. In that case operators 15^*^2 have to 



be computed numerically (for instance using the Method of Fictitious Sources or the 
Differential Method [|T§). 

Around each rod Dn the diffracted part of the field has the following expansion: 
^ F^'^if^^-* {x^n) e*'^'^" where (r„, yj^) are the polar coordinates associated to fiber n. The 

k 

Hankel part of the total field is the sum of the fields diffracted by each rod D„ inside VL: 
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n=l I 



n=l..N 



such that F° 



n=l ^r. 



(5) 



where 



so that there is a hnear operator C 

Conversely, there is a hnear operator 71 = (7^™)^Jf ^ , obtained from the translation 
formula for Bessel functions [T^, such that 7^^ ("^* ~) ~ ^I'n- Operator 71 is injective, i.e. 
it is left invertible, whereas C is surjective, i.e. it is right invertible: C x 7Z = N Id- The 



multi-scattering theory without boundary 
operator 7i such that 



shows that it is possible to define a linear 



7i- 



7ZF' 



(6) 



from which we derive the scattering matrix of the system of rods in the absence of a boundary, 
that is when the medium of permittivity Er fills the entire space: Swb= C7i7Z, which satisfies 
F*^'" = 5^feF*'~. System (§) then rewrites 




Si 1S2 Swb 

Si S2Su]b 




(7) 



from which we derive the expression of the internal and exterior fields from the incident 
field: 

-1 ^-^ 




Swb {id + S2 Swb) Si F*'" 
Si + S2Swb {id + Swb) Si 



(8) 



and the problem is solved. This formulation allows to use distinct numbers Uc (exterior 
cylinder) and nj (fibers) of Fourier-Bessel coefficients. This point is a crucial advantage for 
a low computation time. Indeed, for a given wavelength the interior diffraction problem 
may be correctly described with a small number of Fourier-Bessel coefficients nf whereas 
the exterior problem requires a larger number of coefficient n^. 

In the case of a diffraction problem, F*'+ represents the incident field illuminating the 
structure such as a plane wave, a gaussian beam or a cylindrical wave. The plane wave 
defined by the trihedron (E*, H*, ko) is spatially characterized by its Euler angles {ipo, Oq, Sq) 
where (po, 6q and 5q are respectively called precession, conicity and polarization angle, see 
fig. 2. Therefore, the Cartesian coordinates of kg are given by: 

{/co sin 6*0 cos ipo 
ko sin 6*0 sin cpo 
ko cos 9o 

In this formalism, the case of s (resp. p) polarization is defined by the parameters ^0 = 90 ° 
and 6q = 90° (resp. 60 = 0°). This corresponds to the cases where the unique non zero 
component of the electric field (resp magnetic field) is (resp. Hz). 
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III. NUMERICAL EXAMPLE 



We consider a hexagonal photonic crystal constituted by 19 air holes embedded in a 
dielectric circular cylinder of optical index = 4, see fig. 3. The radii of the holes and 
of the cylinder core are respectively equal to rj = 0.8 and R = 10. The two-dimensional 
structure is illuminated by a plane wave in p-polarization. The transmission coefficient T is 
defined as the fiux of the Poynting vector of the total field collected on a segment situated 
below the cylinder sheath and normalized to the incident energy, see fig 3. The convergence 
of the vector multi-scattering method is studied with respect to the numbers Uf and ric 
of Fourier coefficients respectively used for the series expansion of the interior (fibers) and 
exterior (cylinder sheath) problems. Figure 4 presents the relative error of the transmission 
T versus Uf and ric for an incident plane wave with \/R = 1.70 i.e. in the resonant domain. 
The vector multi-scattering method converges with a relative error less than 0.1% with 
ric > 15 and for nf higher than 3. However, independently of Uc, the number of Fourier 
coefficients must be higher than 4 for a correct convergence of the diffraction problem. 
Therefore, in order to get a good numerical convergence in this wavelength domain, we 
choose ric = 20 and rif = 4. In that case, using a Personnal Computer with a 200 MHz 
processor and with 64 Mo of RAM, the computation time for the transmission coefficient 
for a wavelength and for the structure defined in fig. 3 is about 50 seconds. We have also 
tested our results against that obtained using the fictitious Sources Method [p!5| , p!^ . 

We now study the scattering properties of the PC embedded inside the cylinder sheath 
defined in fig. 3: 

We start with a hexagonal PC constituted by 19 air holes embedded in an infinitely dielectric 
medium of optical index y/e^ = 4 (i.e. without the circular boundary). The PC presents a 
photonic band gap for the interval of wavelengths X/R = [2.0; 2.70] for p-polarization, see 
the bold curve of figure 5. The solid curve represents the transmission T of the PC without 
boundary and doped by a central microcavity (the central air hole is removed) computed for 
p-polarization (6*0 = 90° and 6o = 0°). It appears a resonant wavelength Ar(90°)/-R = 2.13 
inside the photonic band gap associated with a localized mode. 

Let us now compare these results with that corresponding to the structure defined in fig. 3 
(PC with the circular boundary): 

Figure 6 gives the transmission diagram versus the wavelength for a p-polarized incident 
plane wave. The solid curve represents the transmission in the case where the central air 
hole is removed whereas the dashed curve represents the transmission for the perfect crystal 
(with the central hole). This structure also presents a photonic band gap for the interval 
[2.15; 2.7] but slightly shifted toward higher wavelengths in comparison with the PC without 
the circular boundary. We can also remark that though the cylinder sheath is illuminated 
in the resonant domain, the photonic band gap phenomenon provokes the extinction of the 
electromagnetic modes of the exterior cylinder. 

The doped PC embedded inside the cylinder sheath presents a maximum of transmission 
inside the photonic band gap domain for the wavelength A^(90°)/-R = 2.20. The map of 
the modulus of the magnetic field confirms that the resonant wavelength A^(90°)/-R = 2.20 
is associated to a localized mode of the structure, see fig. 7. In conclusion, the resonant 
wavelengths of PCs inside an infinite dielectric medium are shifted toward higher 
wavelengths when PCs are embedded inside a cylinder sheath. These results 
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demonstrate that rigorous numerical computations of finite-size structures are necessary to 
efficiently describe diffraction properties of PC fibers. 

Recent studies have demonstrated that resonant wavelengths associated to localized 
modes strongly depend upon the propagation coefficient 7 (i.e. upon the conicity angle 
60 in the case of scattering problems) PJ20[]. Numerical and theoretical experiments have 



shown that the dependence of the resonant wavelengths with respect of the conicity angle 
is given by: 

A^(^o) = A^(90°)sin^o (9) 

This behavior can be interpreted thanks to a simple model called "model of the infinitely 
conducting cavity" [0. Although the localization of light inside PCs doped by microcavities 
is a global phenomenon due to the photonic band gap effect and the broken symmetry 
of the lattice, both resonant wavelengths and localized modes may be well approximated 
using a local representation of the electromagnetic field. Figure 5 presents the diagram of 
transmission in the case of the doped PC without boundary for 3 distinct conicity angles 90°, 
80°, 70° and for a fixed polarization angle 6q = 0°. The resonant wavelength is shifted toward 
the shorter wavelengths inside the band gap when the conicity angle decreases. In table 1, 
both numerical computations and predicted resonant wavelengths versus the conicity angle 
are compared. The domain of validity of the "infinitely conducting cavity" model depends 
on the strength of the localization of the light with respect to the conicity angle. For example 
in our case, decreasing the conicity angle induced the shift of the localized modes toward 
the lower edge of the band gap. Therefore the localized mode is delocalized out of the 
microcavity hence the model of "the infinitely conducting cavity" is no more valid. This 
behavior explains why equation (|^) must be applied for low inclinations of the wave vector 
with respect to the plane defined by the cross section of the fibers (i.e. for 6 < 60°). Let 
us now study the optical properties of the same PC embedded inside the cylinder sheath 
defined in fig. 3 versus the conicity angle. Figure 8 presents the diagram of transmission 
obtained for 3 values of the conicity angle (90°, 80°, 70°) and for a fixed polarization angle 
5o = 0°. These results lead to two remarks: 

1) When decreasing the conicity angle, a second resonant wavelength appears, 

2) The shift of the resonant wavelength is weaker for a PC fiber than for a PC without 
boundary. 

The apparition of the second resonant wavelength A^(70°)/i? = 2.226 demonstrates that the 
optical interactions of the PC and the cylinder sheath modify strongly the photonic band 
structures. Moreover, the conicity angle dependence of the resonant wavelengths diminishes 
when the PC is embedded inside a cylinder sheath. In that case, the model of the "infinitely 
conducting cavity" cannot be used to compute the shift of the resonant wavelengths for 
grazing incidence. This behavior may be efficiently used for controlling the detunning of the 
localized modes with respect to the inclination of the wave vector. 



IV. CONCLUSION 

We have extended the vector multi-scattering theory of diffraction by parallel cylinders 
to the case of a hard boundary. We have applied this theory in order to study the scattering 
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properties of a PC embedded inside a cylinder sheath. The numerical results have shown 
that PC fibers present complicated photonic band structures with additional localized modes 
for grazing incidence. The exterior circular boundary also attenuates the detunning of 
the resonant wavelengths with respect to the inclination of the wave vector (i.e. with the 
propagation coefficient 7). Matters related to the optical effects of the size of the PC inside 
the cylinder sheath will be discussed in future studies. This theory allows to study the 
transmission properties of 2D PC with inverted contrast but it can also be straightforwardly 
applied to the study of propagation phenomena in photonic crystal fibers, where the modes 
are linked to the resonances of the scattering matrix |jIB|,TP,pT[ Work is also in progress in 
that direction 
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Figures captions: 

Figure 1: Scattering by a set of parallel fibers of arbitrary shape, optical index and 
position. 

Figure 2: Definition of the incident angle (fo, polarization angle 6q, conicity angle in 
the Cartesian system {0,x,y,z). 

Figure 3: Hexagonal PC constituted by 19 air hole embedded in a dielectric circular 
cylinder of optical index ^/e^ = 4. When central air hole (dashed fiber) is removed the PC 
is doped by a microcavity. The radii of the air hole and the cylinder sheath are respectively 
r/ = 0.8 and R = 10. The segment below the structure is used for the computation of the 
transmission coefficient T. 

Figure 4: Relative error of the transmission coefficient T versus n/ and ric- The structure 
is illuminated by a plane wave in p-polarization and for X/R= 1.70. 

Figure 5: Diagram of transmission for the PC defined in fig. 3 but without the circular 
cylinder sheath. The hexagonal PC presents a photonic band gap approximately equal to 
[2.15; 2.70] for p-polarization (bold curve). The solid, dashed and dotted-dashed curves are 
respectively obtained for the incident field parameters: 6q = [90°; 80°; 70°] and for the same 
polarization angle 6q = 0°. 

Figure 6: Logarithm of transmission versus the wavelength for the structure defined in 
fig.3 and for p-polarization: the solid curve is obtained when the PC is doped by a central 
microcavity whereas the dashed curve is computed for the perfect PC. 

Figure 7: Map of the modulus of the magnetic field for the wavelength associated to 
the localized mode A^(90°)/i? = 2.20. 

Figure 8: Diagram of transmission for the structure of fig.3 for 3 conicity angles Oq = 
[90°; 80°; 70°] and for a fixed polarization angle 6o = 0°. 

Table caption: 

Table 1: Comparison between the resonant wavelength versus the conicity angle calcu- 
lated with a direct numerical computation and thanks equation (^. 
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Figure 1 



Conicity angle 0o 


Resonant wavelength 
numericaly computed 


Resonant wavelength 
computated thanks eq.(9) 


90° 


2.12 


2.12 


80° 


2.09 


2.088 


70° 


1.98 


1.992 



Table 1 




V=sin(0o)ez-cos(0o).(cos((po)ex +sin((po)ey ) 
W=-sin((po)ex+cos((po)ey 



Figure 2 



Incident plane wave 
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